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Abstract 

We present the theoretical basis and sample Monte Carlo data for the YFS expo- 
nentiated calculation of e + e~ — > W + W~ — > f\f[ + /2J2 a ^ an d beyond LEP2 energies, 
where the left-handed parts of fi and f[ are the respective upper and lower components 
of an SU2L doublet, i = 1,2. The problem of gauge invariance of the radiation from the 
unstable charged spin 1 W ± is solved in an entirely physical manner. Our formulas are 
illustrated in a proto- typical YFS Monte Carlo event generator YFSWW2, wherein both 
Standard Model and anomalous triple gauge boson couplings are allowed. 
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1 Introduction 



The problem of the precision calculation of the process e + e~ — > W + W~ + 71(7) — > 
Afermions + 71(7) at and beyond LEP2 energies is of considerable interest in connection 
with the verification and tests of the SU2L x U\ model of Glashow, Salam and Wein- 
berg jL] of the electroweak interaction. Indeed, these processes are the primary objective 
of the initial LEP2 physics program, providing as they do both a window on the most 
precise measurement of the W rest mass and a window on the most precise test of the 
fundamental non-Abelian triple and quartet gauge field self-interactions in principle, for 
example. In this paper, we present the theoretical basis of the rigorous Yennie-Frautschi- 
Suura (YFS) @ Monte Carlo approach Q to these processes. For completeness, we shall 
illustrate our results with a proto-typical Monte Carlo event generator which will be exact 
in the infrared regime and will be of leading logarithmic accuracy through 0(a 2 ) in the 
hard radiative regime. A more accurate realization of our results will appear elsewhere ||]. 

Referring now to the results in Refs. || 0, it is clear that there are problems of 
principle in carrying through a manifestly gauge invariant realization of the production 
and decay of massive charged spin 1 boson pairs in e + e~ with the presence of radiation. 
Indeed, some controversy did exist in the literature on just how one should proceed even in 
the stable particle case |5], [|, where for example the current splitting approach of Ref. 
has been questioned as to accuracy and appropriateness in Ref. . We note that recently 
it has been verified |7[] that the approach in Ref. H is indeed accurate enough for the 
requirements of the LEP2 physics program as it is currently envisioned. We will show 
that the YFS theory will afford us an arena in which we can resolve this controversy. 

While we were preparing this manuscript, we became aware of independent and related 
results by Baur and Zeppenfeld (B-Z) || for the problem of quark-pair annihilation into 
followed by a lepton-pair W decay, a process of interest in hadron-hadron collider 
physics for example. We will therefore compare our approach with that of B-Z in what 
follows. We shall see that the two approaches are consistent with one another. Moreover, 
during this same period, we became aware of an independent derivation in Ref. of 
the solution to the gauge invariance problem for radiative corrections to the processes of 
interest to us here. We will therefore use the results of Ref. in what follows, as they are 
in complete agreement with our work insofar as the issue of gauge invariance is concerned. 

In addition, in the original YFS paper ||, there is no explicit discussion of charged 
spin 1 massive radiation. Thus, we will need to extend the the rigorous YFS theory to 
this well. This means that our analysis is of theoretical interest in its own right 

as a study of the infrared limit of massless vector radiation from massive charged vector 
fundamental particles. 

Our work is organized as follows. In the next section, we set our notational conventions. 
In section 3, we derive the extension of YFS theory to spin 1 charged particles. In 
section 4, we derive the corresponding YFS formula for the processes e + e~ — > W + W~ + 
n {l) ~^ 4/ ermions + 71(7) and show that it is manifestly gauge invariant. In section 5, we 
illustrate our YFS formula with Monte Carlo data based on the proto-typical MC event 
generator YFSWW2 flO| , which uses the Born level cross section of Ref. fy] as an input 
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to achieve LL accuracy in the hard radiative regime and exactness in the infrared regime. 
Section 6 contains our summary remarks and, finally, some useful formulae are collected 
in appendices. 



2 Preliminaries 

In this section we review the relevant aspects of our YFS Monte Carlo methods as they 
pertain to the problem of extending our methods to the W + W~ pair production processes 
of interest to us here. In this way we also set our notation and define our kinematics. 




Figure 1: The process e+e — > W + W — ► Afermions = f± + /{ + f% + fy, where \j<), i = 1, 2, 

are SU2L doublets. Here, pa is the 4-momentum of A, A = fi,f-, pt(qx) and ^2(^2) are the 4- 
momenta of e + (e~) and W + (W~) respectively. We use the notation Cl = PlC = |(1 — 75) C for 
any C. 



More precisely, the problem we study herein is illustrated in Fig. Ill, together with the 
respective kinematics: e + + e~ — > W + + W~ +n(j) — > Afermions + ^(7) at CMS energies 
y/s > 2M W , so that we may neglect m 2 e /s compared to one. This corresponds to the 
case of LEP2 and of the NLC (level designs of the NLC now are in progress at several 
l2"fl). Let us focus for the moment on W + W~ pair production part of Fig. 
for the case that the W's are replaced by the fermion pair //, of rest 



laboratories 
In Refs 



13, 14 



masses nif and of charges ±Q/e, we have realized by Monte Carlo methods the process 
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e+ + e 



f + / + 71(7) via the fundamental YFS formula 



dcr = e 



2aSB+2aB 



^ j_ /' r/ {/ ', ; /'_^_ e '." | /'L+7j-/--'/j-^,A-,i-/J 
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where the real infrared function B and the virtual infrared function B are given in Refs. [0, 



13 , 15, 16, 17], and where we note the usual connections 

2a B = 
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for the standard YFS infrared emission factor 

S(k) 



a 



An 2 



"" J ^ pik r/,/,' J 



(2) 



(3) 



if Q f is the electric charge of / in units of the positron charge. Here, the ". . . " represent 
the remaining terms in S(k) obtained from the one given by respective substitutions of 
Qfj Pij Qf'i Qi with corresponding values for the other pairs of the respective external 
charged legs according to the YFS prescription in Ref. @ (wherein due attention is taken 
to obtain the correct relative sign of each of the terms in S(k) according to this latter 
prescription) and in Ref. |13|, [14|], f ^ e, f — f. 
The YFS hard photon residuals /3j in (1) 



i — 0, 1,2, are given in Refs. |T3|, 13 for 
YFS2, YFS3 so that these latter event generators calculate the YFS exponentiated exact 
0(a 2 ) LL cross section for e + e~ — > // + n(j) with multiple initial, (initial+fmal) state 
radiation respectively using a corresponding Monte Carlo realization of (1). In the next 
sections, we use explicit Feynman diagrammatic methods to extend the realization of (1) 
in YFS2, YFS3 to the corresponding Monte Carlo realization of the respective application 
of (1) to e + e~ — > W + W~ + 71(7) — > Afermions + 71(7). 



3 YFS theory for massive charged vector particles 

In this section, we present the required formulas for extending the exact O (a 2 ) LL Monte 
Carlo realization of the hard photon residuals /3 n in YFS2, YFS3 |L3], Q to the corre- 
sponding Monte Carlo realization for the W /+ W / ~-pair production processes. We begin by 
deriving the respective YFS real and virtual infrared functions for these latter processes. 

Referring to the kinematics in Fig. [I] and to the definition of the YFS infrared functions 
S(k), S(k) in Ref. 0, we see that to extend the infrared YFS calculus to W's it is enough 
to show that in the respective infrared regime for an emitted photon of 4-momentum k the 
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amplitude for spin 1, charge Qw& and mass My/ for the emitting particle is related to that 
for spin 1/2, charge Qfe and rest mass rrif by the substitutions of the respective Lorentz 
group representation factor, the charge, the corresponding 4-vectors and radiationless 
(Born) amplitude. From this result, it is immediate that the formulas given in Ref. JJ] 
for for S(k), S(k) for spin 1/2 hold also for spin 1 with the corresponding substitution of 
charges and massive 4-vectors. Let us now establish this correspondence of the infrared 
regimes of spin 1/2 and spin 1 massive, charged particles. 





Figure 2: Infrared emission of a photon of 4-momentum k by (a) an incoming fermion of mass rrif, 
charge Qfe and 4-momentum p, (a') an incoming spin 1 charged boson of mass My/, charge Qy/e 
and 4-momentum p, (b) an incoming anti-fermion of rest mass rrif, charge —Qfe and 4-momentum 
p, and (b') an incoming spin 1 charged boson of rest mass My/, charge — Qy/e and 4-momentum 
p. Here, the components of p of course must change in passing from the fermion cases to the vector 
boson cases if the incoming lines are on-shell, for example. 

The relevant situations are illustrated in Fig. |2], where in (a) we have the emission of 
a photon of 4-momentum k by an incoming fermion of charge Qfe, rest mass rrif, and 
4-momentum p to be compared with the analogous situation in (a') where we have the 
emission of a photon of 4-momentum k by an incoming spin 1 vector boson W~ of charge 
Qw&i rest mass My/ and 4-momentum p. In (b), we have the emission of a photon of 
4-momentum k by an incoming spin 1/2 anti-fermion of charge —Qfe, rest mass rrif, and 
4-momentum p to be compared with (b') wherein we have the emission of a photon of 
4-momentum k by an incoming spin 1 vector boson W + of charge — Qy/ e i res t mass My/ 
and 4-momentum p as well. From the standard Feynman methods, for Fig. 0(a), we have 



4 



the amplitude, for / = e for definiteness, 



M 2a = ■■■ -r, % —— (-iQ e erf)u(p) 

f — m e + ie 

iQ,<) {i>- jt + m e h»u(p) 



-2kp + k 2 + ie 
(Q e e) 



k 2 — 2kp + ie 



[2p» -k» + ia» a k a ]u(p) 



^ = - f- ( y +it |y - t ' w ' 



where we define 



lim A(k) = \imA(k), 

IR k~>0 

for any function A(k). For Fig. 0(a'), we get the corresponding result 
M2a ' ~ (p* -M 2 W + ie) {%Qwe) 



[g a ,p(2p - ky + gl,{-p + 2k)„ + g<t(-p' - 2k) a ,]e p _{p) 



_ (Qwe)(g^-p' a p' a /M 2 w ) 
~ (k 2 - 2k P + ie) ^(2 P -k) 

+ 2kpgZ,-2k a ,g»-p' a ,g%}et(p) 

Similarly, for Fig. 0(b), we get the infrared limit as 

M 2b = v(p){-iQ e e)Y — — — • ■ ■ 

— p+ ft — m e + ie 

v(p)Q e e^(- ft + m e ) 

k 2 — 2kp + ie 

_ Qeev(p)[-^ p 1 - fa* + (p' + m e )Y + 1(^ + 7'* #) + K7" J£- #7")] 

k 2 — 2kp + ie 
Q e ev{p)[-2p^ + k»- ia^ka] 



limA^b 



k 2 — 2kp + ie 
v{p)Q e e{-2p^ + 



ir " /c 2 - 2kp + ie 



whereas for Fig. §(b') we compute the infrared limit as 

_ H)(9^-P' a "p'"/M^ 

Mw ~ (p' 2 -M 2 W + ie) {lQwe)[ 

gp a (-2p + kf + g£(p' + 2% + $(p - 2k) a }e a + (p) 

{Qwe){g a "^v' a "v' P /M^) 



(k 2 - 2kp + ie) 
+ 2kpgZ-2k a g% + p' p gZ}el(p) 



[gp<x(-2p + k) 



This shows that the stated correspondence holds. 

It follows that we obtain the YFS infrared functions |2|] S(k) and S(k) for real and 
virtual soft photon emission from W + '~-\ines by substituting the respective mass M w into 
the corresponding expressions for these functions for emission from e +,_ -lines: 



(8) 



where the subscripts indicate the respective YFS infrared functions for ee, ee, ee, W~W + , 
W~W~ and W /+ VI^ + in the obvious manner. 

One important point needs to be discussed before we turn to the application of (8) 
to W ± pair production. This concerns the so-called Coulomb effect [pl| — the enhanced 
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1/(3 behavior of the 0(a) virtual correction to the Born cross section due to the exchange 
of k — > virtual photons, where j3 is the CMS velocity of one of the TV's. Since the 
YFS virtual infrared function 2aB\y-w+ = J [d^kjik 2 — m 2 + ie))S(k)w-w+, where m 7 is 
our photon infrared regulator mass, describes precisely this regime of the virtual photon 
phase space as well, we need to remove this Coulomb effect from B w - W + so that it can be 
treated via the methods of Ref. ]18[ as accurately as one desires without double counting 
it. This we do by defining the analytic subtraction 

B' W . W+ (P) = B w - w+ {(3) - diPt ~ P) hmPB w - w+ (P) 

= B w - w +((3) - jpO(Pt ~ 13), (9) 

where 9 is the usual step function and here we determine (3t, the transition velocity which 
separates the non-relativistic regime from the relativistic one insofar as the Coulomb cor- 
rections are concerned, by the requirement that the 0(a 2 ) Coulomb correction, (^) 2 /12, 
is less than 0.03% for (3 > (3 t — this gives (3 t = 0.382. For (3 > (3 t , the Coulomb cor- 
rection series is a well-behaved part of the usual perturbation series and does not need 
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special treatment in our analysis. Our requirement that the 0(a 2 ) Coulomb correction 
stay below 0.03% for (3 > f3 t is determined with an eye toward an ultimate goal of 0.1% 
total precision on our calculations — this goal is however beyond the scope of the current 
paper. It is apparent that B' w _ w+ and B w -w+ contain the same infrared divergences so 
that we may use B' w _ w+ in our YFS exponentiation algebra without any change in the 
cancellation of infrared singularities to all orders in a so that such a use of B' w _ w+ is 
fully justified by the original YFS arguments - the resulting cross sections are fully 
independent of the substitution of B' w _ w+ for B w - W + when the theory is summed to all 
orders in a as it is proved in the original YFS paper. Thus, we will make this substitution 
here and treat the Coulomb effect entirely according to the methods in Ref. fl8fl . It can 
therefore be seen that our procedure for arriving at a smoothe transition, in our complete 
cross section, within the limits of its physical precision tag, between the Coulomb dom- 
inated j3 — > regime and the relatavistic (3 — > 1 regime of the W + ~ charge form factor 



is entirely consistent with the smoothe interpolation of Schwinger in Ref. between 
the analogous regimes of the respective charge form factors for spin and spin ~ massive 
charged particles. 

We now turn to the application of the results in this section to the realization of our 
YFS Monte Carlo approach to multiple photon radiative effects in e + e~ — > W + W~ + n( / y) 
at LEP2 and NLC energies. 



4 Gauge invariant YFS Monte Carlo for e + e —> W + W + 



In this section we apply the results of the preceding section to develop a Monte Carlo 
event generator, YFSWW2 [ID"], which realizes the YFS exponentiated multiple photon 



radiation in the process e + e~ — > W + W~ + ^(7), where we will also allow the W's to 
decay to four-fermion final states. In the development presented here, we shall work to 
the leading log /3 level in the respective YFS hard photon residuals /3 n in (1) as it is 
applied to VF-pair production via the substitutions in (8). 

Specifically, on using the results in (8), we arrive at the representation, for the process 
e + e~ — > W + W~ + 72(7) — > fi + f[ + f 2 + J 2 + n(j), of the fundamental YFS cross section 
formula 



A n — p 2aReB+2aB J_ / TT ^± f 

^n\J 11 ife? J (2vr) 4 

n=0 J j=l J J v ' 



iy(pi+qi-P2~(i2-Y.j kj)+D 



dPpfrdPpji d 3 p f d 3 p j 2 
p n (kx, . . . , k n ) q t. 7; 7; , (10) 
P° fl P%P%P% 

where, referring to the kinematics in Fig. [I], we have the identifications p<i = pf 1 +Pf, 92 = 
Pf +Pf 2 , for the W + , W~ 4-momenta, respectively. The YFS infrared functions B(K max ), 
D, and B, by (8), are obtained from the the corresponding ones for the process e + e~ — > 
f + f + n(j) as given in (§) and in Refs. [|l|, [I7j via the substitutions: (Qf,pf,rrif) — > 
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{Qw-,Q2,Mw), (Qf,Pf, m f) ~ > (Qw+,P2, Mw)- Their analytical representations are given 
in Appendix A. The hard photon residuals f3 n now contain both the production and the 
decay of the Ws, which may of course occur either on or off the W mass-shell. We will 
work to the /3q level and we have the identification 

do Bom 



(11) 



where c?f2 + (_) is the differential decay solid angle of fiify) i n the rest frame for 

example. Here, we take the respective Born cross section, da Bom, from Ref. [II] for 
definiteness. The result ([TU]) has been realized via the YFS MC methods of two of us, see 
for example Ref. 



PJ, in the program YFSWW2 [1C . 



-Pi 



qi 




+ 



Figure 3: The imaginary parts necessary to include in the amplitude for e + e~ — » W + W~ — > 
Afermions = f\ + /{ + f 2 + f 2 in order to maintain gauge invariance in the presence of radiation 
by the Ws themselves. Here, the notation is identical to that in Fig. |l| and the vertical dotted line 
indicates the standard Bjorken-Landau-Cutkosky isolation of the respective imaginary part. 

Before we illustrate the type of predictions which we make with YFSWW2 for the 
LEP2/NLC physics scenarios, let us address one further important theoretical point con- 
cerning the fundamental result ([10|) . Specifically, the reader may note that the W ± decay 
width Ty/ does not appear in the YFS infrared functions B and B as they are given via 
the result (8). Yet, evidently, when a W ± of 4-momentum p' = p + k emits radiation of 
4-momentum k, the propagator denominator that would most naively be present imme- 
diately preceding the emission vertex would be D\y(p') = p' 2 — M"^ + ip' 2 T\y / 'Mw and for 
p 2 = M^r this does not reduce to the YFS infrared algebraic form k 2 + 2kp. The atten- 
dant problems with electromagnetic gauge invariance are now well-known || |6|, || |9| . The 
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solution to this apparent dilemma is already presented in Refs. [[| ||. Whenever the W ± 
radiate, one can include their decay width in their propagation provided one also allows 
the radiated photon to couple to charges in the graphs which generate the non-zero width 
to its respective order as well, see Refs. || || for a detailed discussion of these graphs, 
which are illustrated in Fig. |^ for the case that T w is computed to 0(a) for example. The 
net effect, as explained in detail in Refs. |!| and as we have independently verified, is 
that the emission vertex for the photon in question is multiplied by a factor which replaces 
the pre-emission denominator Dw(p') = p' 2 — + ip' 2 Tw/Mw for p 2 = with the 
desired YFS infrared algebraic form k 2 + 2kp, thereby maintaining electromagnetic gauge 
invar iance. 

We turn now to sample Monte Carlo data from YFSWW2 in the LEP2 and NLC type 
energy regimes. This we do in the next section. 

5 Sample Monte Carlo data for e + e~ — > W + W~-\-n(j) — * 
Afermions + 71(7) 

In this section, we present sample Monte Carlo data for the process e + e~ — > W + W~ + 
n (l) 4:fermions + ^(7) in the LEP2/NLC type energy regimes. We have used our 
Monte Carlo event generator YFSWW2 as presented above and we have had in mind in 
particular the illustration, among other things, of the effect of the contribution to 
the YFS real and virtual infrared functions in (|10D, as this effect has not been presented 
elsewhere. 



Ecu [GeV] 


ISR 


ISR ± Coul. corr. 


ISR ± Coul. corr. ± Y'-corr. 


175 


0.4906 ± 0.0002 
0.4898 ± 0.0002 


0.5046 ± 0.0002 
0.5037 ±0.0002 


0.5053 ±0.0002 
0.5048 ±0.0002 


190 


0.6060 ±0.0007 
0.6034 ±0.0007 


0.6193 ±0.0007 
0.6166 ±0.0007 


0.6217 ±0.0007 
0.6195 ±0.0007 


205 


0.6359 ± 0.0008 
0.6315 ±0.0008 


0.6480 ± 0.0008 
0.6436 ± 0.0008 


0.6514 ±0.0008 
0.6476 ± 0.0008 


500 


0.2910 ±0.0003 
0.3538 ± 0.0004 


0.2946 ± 0.0003 
0.3582 ± 0.0004 


0.2970 ±0.0004 
0.3591 ±0.0004 



Table 1: The results of the 10 5 (except for Ecm = 175 GeV, where it is 10 6 ) statistics sample 
from YFSWW2 for the total cross section a [pb\. The upper results at each value of energy are for 
Standard Model couplings constants, while the lower ones are for anomalous couplings constants 
(5k = 5X = 0.1). See the text for more details. 



Specifically, we will use the 4-fermion final state f% — c, f[ — s, f 2 



f'2 = e for def- 



initeness. With an eye toward studies of the triple gauge boson couplings in the standard 
model SU2 x U\ theory of Glashow, Salam and Weinberg [|I|, we further follow the notation 
of Ref. |Tl| and compute our results for no anomalous 3-gauge boson couplings and for 
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the values 5 k = 0.1, 5X = 0.1 for the deviations of the coupling parameters k, X away from 
their tree-level standard model expectations of 1.0, 0.0, respectively, in the notation of 



Ref. 0. This we do for three LEP2 type CMS energies y/a = 175, 190, 205 GeV and for 
the NLC-type energy 500 GeV. The cross sections which we obtain from our simulation 
with YFSWW2 were each determined from 10 5 events (except for 175 GeV, where we had 
10 6 events) and for the full solid angle acceptance for each of the 4-fermions in the final 
state. The results were computed for both the pure initial state YFS exponentiated /So- 
level case, as in Ref. ]20| , referred to as the ISR case in the following, and in the case when 
the W ± contribution to the soft YFS exponentiated radiative effects is treated exactly, 
which we refer to as the ISR+Y' case in the following. Furthermore, for definiteness, and 



clarity, the effect of the Coulomb correction, as given in Ref. [18], is also illustrated in our 
results, both for the simple YFS form factor ISR case and for the full YFS form factor 
ISR+Y' case — the presence of the Coulomb correction is indicated by "Coul." in the 
following. The corresponding results are illustrated in Table [1]. 

What we see is that the effect of the full YFS form factor is at the level of 0.14%, 
0.39% and 0.52%, respectively beyond the usual Coulomb effect at the LEP2 energies 
175 GeV, 190 GeV and 205 GeV and is at the level of 0.81% beyond the usual Coulomb 
effect at NLC energies for the case of the SM couplings, giving a total effect beyond 
initial state radiation of 3.00%, 2.59%, 2.44% and 2.06% for the CMS energies 175 GeV, 
190 GeV , 205 GeV and 500 GeV , respectively. For the case of the anomalous couplings of 
5k = 5X = 0.1, the corresponding results are 0.22%, 0.47%, 0.62% and 0.25%, respectively 
beyond the usual Coulomb effect for total corrections of 3.06%, 2.67%, 2.55% and 1.50%, 
respectively for the same CMS energies. Thus, we see that at LEP2 and NLC energies, 
the full form factor does indeed modulate the effect of the anomalous couplings; this 
enhances its importance at both LEP2 and NLC energies. Since the targeted accuracy 
of the theoretical precision for the LEP2 lyi^-pair production cross section is 0.5% 0], 
evidently, the full form factor effect is very important for LEP2 physics scenarios. 

Indeed, we have looked into the manifestation of these effects in the production 
angular distributions. We show in Figs. |4], [| four respective differential distributions for 
the total cross sections in Table [l] for the LEP2 energy of y/s = 190 GeV and for the NLC 
energy of \/s = 500 GeV, where we feature the W ~ production angle distribution in the 
CM system for all three SM coupling scenarios in Table |l] and with the final "ISR+Coul. 
corr.+Y'-corr." scenario for the anomalous coupling case. We see in these figures that 
the full formfactor effect modulates the distributions most strongly near their peaks, near 
the forward direction for the W~ when the incoming electron direction is used as the 
reference direction. As the anomalous couplings modulate these distributions over a large 
range of the respective production angles, we see that the smaller (larger) values of the 
full formfactor effects in Table [1] for the anomalous cases relative to the SM cases for 
NLC (LEP2) energies is consistent with the shapes of the distributions in Figs. |4], |5]. Also 
evident in the figures is the more pronounced anomalous coupling effect at NLC energies, 
as expected. 

We end this section by noting that we have also implemented the full YFS form factor 
effect and the anomalous couplings as well in an unpublished version of the program 
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da r pb 



da r 



0.5 



0.4 - 



0.3 - 



0.2 - 



0.1 '-{ 



0.0 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



— ISR 

x ISR + Coul. corr. 

* ISR + Coul. corr. + Y'-corr. 

o ISR + Coul. corr. + Y'-corr. + An. c. c. 




i i i i i 



0.5 1.0 1.5 2.0 2.5 3.0 

6 [rad] 



Figure 4: W angular distributions for Ecm = 190 GeV . See the text for more details. 



KORALW [[2(J of three of us (M.S., S.J. and W.P.) and we have checked that the results 
from YFSWW2 and from this new version of KORALW are in agreement within the 
statistical errors of the simulations. This is an important cross-check on the results in 



this paper it will be presented in detail elsewhere [21 . 



6 Conclusions 

In this paper we have developed the YFS theory of charged, spin 1 bosons in the presence 
of possible non-zero widths. We have applied our theory to the LEP2/NLC process 
e + e~ — > W + W~ + n( / ~f), allowing for the W-pair to decay to 4 fermions. The result is the 
description, via the Monte Carlo event generator YFSWW2, of the respective n(j) effects 
on an event-by-event basis, in which the infrared singularities are cancelled to all orders 
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Figure 5: W angular distributions for Eqm = 500 GeV . See the text for more details. 



in a. 

Specifically, in our realization of the YFS theory for charged spin 1 bosons, we have 
maintained the electromagnetic gauge invariance of the SU2L x U\ theory following the 
works in Refs. || ||. In addition, we have also avoided any doubling counting of the so- 
called Coulomb effect by removing it analytically from the YFS virtual infrared function 
B. This resulted in the definition of a new YFS virtual infrared function B' . We have 
illustrated our calculations with explicit Monte Carlo data at the LEP2 CMS energies 
v/i = 175 GeV, 190 GeV, and 205 GeV and at the NLC energy = 500 GeV, for both 
Standard Model and anomalous VWW vertices, V = 7, Z. We find in all cases that the 
effect of the radiation by the W ± themselves is important, both in the production angular 
distributions and in the over-all normalization. In our YFS Monte Carlo realization of 
this effect, we have worked to the leading (3 level in the current analysis. The higher 
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order hard photon residuals j3 n , n > 1 will be considered elsewhere [Q in this connection. 

In summary, our Monte Carlo event generator YFSWW2 now calculates, on an event- 
by-event basis, the multiple photon effects in the process e + e~ — > W + W~~ + 71(7) — > 
Afermions + 71(7) and includes, for the first time ever, the effects of the radiation by 
the themselves in the respective YFS exponentiated soft photons, without doubling 
counting the so-called Coulomb effect and without spoiling the electromagnetic gauge in- 
variance of the SU^l x U\ theory. This program is available from the authors at the WWW 
URL |http:/ /enigma.phys. utk.edu/pub/YFSWW /| and we look forward with excitement 



to its application to imminent LEP2 data. 
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A YFS infrared functions 

In this section we present some analytical representations of the YFS infrared (IR) func- 
tions corresponding to emission of virtual and real photons for the VF + W _ -pair production 
in the e + e~ annihilation. An important feature of these representations is that they are 
stable and fast in numerical evaluation. Thus, they are particularly suited for Monte 
Carlo implementations. 

A.l Virtual photon IR function for s-channel 

The s-channel virtual photon YFS IR function dtB(s) for any two charged particles with 
four-momenta pi,p 2 and masses m 1; m2 reads 



a 



1, Ml + P) .\, m 7 , W,K1 + P) , rnj-m 2 2 mi 



2a^B(s, m u m 2 ) = -{ - In — — - 1 In -2- + £Z1 i n LI + _^ £ ] n _± 

7r \p mitn2 J m\m2 s m 1 ro 2 2s m 2 



1 

+ - 
P 



7r In k — In t, In 

2 m\ m\ 2 wij + p(l + p) 

, • , 2 W \ t • ( 2 PP 



m\ + + p) J \m| + + p) 

where 

P = PiP2, s = 2p + m\ + m 2 2) 



(12) 



'-V 1 -^)' (13) 

and m 7 is a fictitious photon mass used to regularize the IR singularity. 
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A. 2 Virtual photon IR function for t and u -channels 

The i-channel virtual photon YFS IR function dtB(t) for two charged particles with four- 
momenta p\,P2 and masses m, M, where m <C M, \t\, reads 

a J / _ \t\ : K A ^ m 2 | ( f u \t\ 



2* mt , m , M) = - ^ta-LL, + ln C - Ij In J + | ^.n -LL, + >n C 

1 , 1*1 , 1*1 , M A 1*1 , > C - 3 

- in — - in — - — in — m h in C H 

2 m 2 M 2 m \ mM 2 

h f( to sl? + 5 b,f ) + "*S)" ' (14) 



where 

C = l + ^p * = (pi-p 2 ) 2 - (15) 

The w-channel IR function dtB(u) can be obtained simply by replacing t — > u in the above 
formula. 

A. 3 Real photon IR function 

The YFS IR function B corresponding to the emission of real photons with energy E 1 < 
K max in a process involving any two charged particles with four-momenta p±, P2 and masses 
777.1,777.2 can be expressed as 



a A fx(l + p) \, 4^ a2 , 1 , 1+A , 1 , l + # 



2aB( Pl ,p 2 - 1 K max ) = -{ -ln " v r/ - 1 In — + — In — ^ + — In 



2 



7r I \p mim 2 / m 2 2/?i 1 — f3\ 2(3 2 1 — A 



2 



+ M 4 (pi,p 2 )>, (16) 



where $ = yl — m 2 /E 2 , and /i and p are defined in Eq. (13). The most complicated part 
of the above expression is the function Ai(pi,p 2 )- It can be expressed as a combination of 
some number of logarithms and dilogarithmsQ (see also Ref. |22| for a similar calculation) 



y/(Q 2 + oo 2 )(Q 2 + 5 2 ) L VA 2 + Q 2 + A 

+ (17) 



using some identities we managed to reduce the number of dilogarithms to 8 only! 
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where 



and 



(v -yk)(v -yi) 



1 



Y( V ) = z u ( v ) + z 21 ( v ) + z 32 ( v ) - z u ( v ) + -xlMxlt^), 



Zij{ri) = mu 2 



Vj - Vi 



1 



In 2 



Vo = 

1/1,2 
1/3,4 



m 



2' 



»7i 



E\ 



7 A 2 + Q 2 , 



a/a 2 + Q 2 ' 1 



A 



v / A T TQ 2 + fi + 



^± y/(g 2 + ^ 2 )(g 2 + (5 2 ) 

v'A 2 • O 2 A 



where we used the following notation 



A — E\ — E 2 , Q — Ei -\- E 2 , 
5 = mi — m 2 , w = mi + m 2 , 

12 /„ „ \2 



Q 2 = -(pi-P2Y 



(18) 



(19) 



(20) 



The only approximation used in deriving the above formulae is m 7 <C K max . 

As one can easily check the dependence of the above functions on the IR regulator m 7 
cancels out in the sum 2a$lB + 2aB which is used to construct the YFS formfactor. 
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